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The interplay of symmetry and topological order leads to a variety of distinct phases of matter, the Symmetry 
Enriched Topological (SET) phases. Here we discuss physical observables that distinguish different SETs in the 
context of Z 2 quantum spin liquids with SU(2) spin rotation invariance. We focus on the cylinder geometry, and 
show that ground state quantum numbers for different topological sectors are robust invariants which can be used 
to identify the SET phase. More generally these invariants are related to ID symmetry protected topological 
phases when viewing the cylinder geometry as a ID spin chain. In particular we show that the Kagome spin 
liquid SET can be determined by measurements on one ground state, by wrapping the Kagome in a few different 
ways on the cylinder. In addition to guiding numerical studies, this approach provides a transparent way to 
connect bosonic and fermionic mean field theories of spin liquids. When fusing quasiparticles, it correctly 
predicts nontrivial phase factors for combining their space group quantum numbers. 
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I. INTRODUCTION 


In contrast to conventional phases that are distinguished by 
Landau order parameters, topologically ordered states with 
emergent anyonic excitations remain distinct even in the ab¬ 
sence of symmetry. However, the presence of symmetries, 
which is natural in most physical contexts, leads to further dis¬ 
tinctions, the so called symmetry enriched topological phases 
(SETs). Well known manifestations include fractional charge 
of anyons in the fractional Quantum Hall states (and in frac¬ 
tional Chern insulators) and spin-charge separation in quan¬ 
tum spin liquids. Recently, rapid progress in the theoretical 
understanding of SETs is being made nun. This is partly 










2 


driven by conceptual advances in the related theory of strongly 
interacting Symmetry Protected Topological phases (SPTs), 
where despite the absence of anyon excitations in the bulk, 
nontrivial edge excitations emerge SEMI. What would 
be particularly welcome at this stage, is a physically well mo¬ 
tivated example of an SET. To make progress in this direction 
we will need to understand how different SET phases can be 
distinguished. 

Further motivation for studying measurable characteristics 
of SETs is the recent progress in the search for quantum spin 
liquids in frustrated magnets, both in experiments and in nu¬ 
merics. A number of S= 1/2 quantum magnets in 2D and 3D 
frustrated lattice have been identified, which appear to evade 
magnetic order 1 E 9 I [20tl . These include the S= 1/2 Kagome 
material, Herbertsmithite, which shown no sign of ordering 
down to temperatures that are a thousand times smaller than 
the exchange constant [21]. While clear cut evidence of an en¬ 
ergy gap in these materials is yet to emerge, a requirement to 
be considered topologically ordered, this may be an extrinsic 
effect due to impurities (although other explanations have also 
been suggested li22U24ln . Meanwhile, extensive density ma¬ 
trix renormalization group (DMRG) simulations of the nearest 
neighbor Kagome Heisenberg antiferromagnet (KHA) indi¬ 
cate (i) a gapped ground state that respect all symmetries [25j] 
and (ii) a topological entanglement ll26l 12711 entropy of log 2. 
In reference ||28il it was argued that the ground state must pos¬ 
sess Z 2 (toric code) topological order to be compatible with (i) 
and (ii). An important open question is the identification of the 
precise phase of matter, i.e. the SET, realized by the Kagome 
antiferromagnet. While a complete solution would necessitate 
extensive numerical input, and is beyond the scope of this pa¬ 
per, here we will relate SET phases to physical properties that 
are readily measurable in numerical simulations. 

A prerequisite is a measurement of the topological order it¬ 
self. Entanglement entropy j29l [301 provides one signature, 
although it does not uniquely specify the topological order. 
A complete characterization is obtained, from either the en¬ 
tanglement spectrum |3li| ( in certain cases), or the S and T 
matrices D2l-l38ll some of which are well suited to numerical 
calculations 1391 - 14211 . Related techniques can be used to diag¬ 
nose 2D SPT phases l43ll . 

Here we discuss measurable properties that distinguish 
SETs. For the reasons above we focus on the S= 1/2 KHA, 
assuming Z 2 (toric code) topological order, which has a pair 
of emergent S= 1/2 spinons (one bosonic and one fermionic) 
and a vortex (vison). The different SETs differ in their re¬ 
alization of space group symmetries and their interplay with 
time reversal. 

We consider systems on both finite and infinite cylinder ge¬ 
ometries, which are well suited to DMRG calculations. We 
show that 1) the many-body quantum numbers of a finite- 
cylinder ground state under space group symmetries such as 
reflection, translation etc. provide a powerful diagnostic of the 
underlying SET and 2) when viewing an infinitely long cylin¬ 
der as a ID spin chain, the ID SPT order of the spin chain for 
various geometries and topological sectors completely deter¬ 
mines the SET order, at least within the space of mean-field 
parton ansatz. 


A reflection quantum number is a probe of quantum entan¬ 
glement: if a state is odd under reflection, then the two halves 
of the system are entangled. But how can quantum numbers 
of a global symmetry, which are always integral, be related to 
symmetry fractionalization? As a simple example - consider 
creating a pair of identical anyons from the vacuum. Say that 
the pair is related to one another by a symmetry (such as a 
reflection or rotation). If the excitations carry charge, a unit 
charge for this state actually implies half charge for each ex¬ 
citation since they are constrained by symmetry. This is one 
of the key ideas that we will exploit. Its implementation is 
more involved when the second symmetry is not charge, but 
also a space group (or time reversal) symmetry. Neverthe¬ 
less such arguments establish the relative quantum numbers 
between different topological sectors. 

To relate our results to an established classification scheme, 
we use a specific model of SETs obtained by a parton decom¬ 
position of the spin operator into bosonic (Schwinger bosons) 
or fermionic (Abrikosov fermions) partons. Symmetry frac¬ 
tionalization is encoded in the Projective Symmetry Group 
(PSG) fE m, mh, which determines how the partons trans¬ 
form under symmetry. A parton mean field theory combined 
with projection leads to a spin wave-function whose quan¬ 
tum numbers in each topological sector are completely de¬ 
termined, which reflect the underlying SET. In particular we 
show that the Kagome spin liquid SET can in principle be 
uniquely determined by measurements on one ground state in 
a few different finite-cylinder geometries. This information 
is numerically superior to the relative quantum numbers be¬ 
tween topological sectors, since DMRG numerics on KHA do 
not obtain all topological sectors in the finite sized systems 
studied. 

A different perspective on our approach is to view it as a 
‘dimensional reduction’ in which we view a 2D SET phase 
on a cylinder as a ID SPT. The nature of the ID SPT de¬ 
pends on the topological sector being studied, the geometry 
of the cylinder, and the SET. In addition to its utility as a di¬ 
agnostic in numerics, our approach provides a theoretical tool 
to study connections between different representations of the 
same SET. The ID SPT invariants for the degenerate ground 
states (which are labeled by quasiparticles of the topological 
order) are shown to follow the same multiplicative law as the 
fusion rules for the Abelian quasiparticles. This allows us to 
determine from the PSGs of two anyons types in the Z 2 topo¬ 
logical order the PSG of the third anyon type, which is found 
to obtain nontrivial phases in certain cases. Some of these 
were not previously known H and serve to correctly relate 
bosonic and fermionic mean field states on the Kagome lat¬ 
tice (H]. In particular this leads us to equate two popular states, 
the Q 1 = Q2 Schwinger boson state ll46h and the Z 2 [0,7r]/3 
fermionic mean field state [5, 47,4$]. 

In addition, we show how dimensional reduction can be 
used to completely identify the four topological sectors of a 
cylinder; this is highly useful for DMRG studies, and does 
not require simultaneous knowledge of all four ground states. 
In particular, we have found a ID SPT invariant that distin¬ 
guishes between the bosonic and fermionic spinon. 

Earlier work employed a similar dimensional reduction ap- 
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proach in the case of internal symmetries with projective rep¬ 
resentations, |4g|] and here we find it to be much more gen¬ 
erally applicable in the presence of space group symmetries. 
Related work specializing to the case of just translation sym¬ 
metry has recently appeared fsoll . The connection between 
ID SPT invariants and global quantum numbers was previ¬ 
ously noted mm, and other works have utilized global 
many body quantum numbers to identify topological phases 

mm. 


H. REVIEW OF Z 2 SPIN-LIQUIDS 

According to the arguments of Hastings, Oshikawa, Lieb, 
Shultz and Mattis, Il54l - l56ll a quantum magnet with half¬ 
integer spin per unit cell is either gapless, breaks a symme¬ 
try, or is a gapped ‘spin - liquid’ with emergent anyonic ex¬ 
citations. In the latter case, the simplest possibility consistent 
with time-reversal is the Z 2 ‘toric code’-type spin-liquid. [ Hi 
This phase has emergent 5=1/2 excitations, the ‘spinons,’ 
even though a truly local excitation (the magnons) must carry 
5=1. These emergent spinons are anyonic excitations with 
non-trivial braiding and statistics. The Z 2 spin-liquid has four 
anyon types: the local excitations (‘1’); the bosonic spinon 
(‘6’), which carries 5 = 1/2; the vison (V), which be¬ 
haves like a 7r-flux for the spinon; and the fermionic spinon 
(‘/’), formed from the composite of / = bv, and which 
also carries 5 = 1/2. Each particle is its own anti-particle, 
v 2 = f 2 = b 2 = 1 (hence ‘Z 2 ’). The braiding and statis¬ 
tics of the Z 2 spin-liquid are equivalent to Z 2 gauge theory 
(the ‘toric code’). In the language of Z 2 gauge theory, b is the 
electric charge e; v is the magnetic flux m; and / is the dyon 
/ = em composed of flux and charge. 

The simplest SET aspect of the Z 2 spin-liquid is its behav¬ 
ior under SO(3) spin rotations: 1 and v carry integer spin, 
while the spinons 6, / carry half-integer spin. The half-integer 
spin carried by the spinons is ‘fractionalized’ because any lo¬ 
cal excitation of the constituent 5 = 1/2 spins transforms as 
5 = 1. We can look for additional SET distinctions based 
on the transformation properties of anyons under space-group 
symmetries, which is the subject of this work. 


A. Minimally entangled states 

The Z 2 spin liquid has a 4-fold topological ground state de¬ 
generacy on the torus or cylinder. Throughout this work we 
rely on a special basis for the ground-state manifold called 
the ‘minimally entangled states’ (MES).|33]] To construct the 
MES basis, let x run along the infinite length of the cylin¬ 
der, and let denote the adiabatic process in which a pair 
of anyons a/a are created from the vacuum and dragged in 
opposite directions ±x out to infinity. The process T" re¬ 
turns the system to the ground state, so is a unitary operation 
in the ground-state manifold. We say that ‘threads any¬ 
onic flux a’ through the cylinder. Since the model is Abelian, 
F". F x ex F"' ! \ where a ■ b denotes the fusion of Abelian 
anyons. 


The MES basis is the unique basis in which J 7 " is realized 
as a permutation of the basis states for all a £ {1 ,b,v, /}. F x 
acts as a permutation in the MES basis because each MES has 
definite topological flux threading the cylinder (its ‘topologi¬ 
cal sector’). 

In previous discussions of the MES it is often assumed 
that each of the four MES can be uniquely identified with an 
anyon type l,b,v, /, so that each MES can be labeled with 
an anyon type |a). Then threading topological flux is realized 
as |a) = | b ■ a). This is the case for Z 2 spin-liquids on 
even circumference cylinders, but for an 5 = 1/2 model on 
an odd circumference cylinder there is a subtlety that arises 
because the MES double the unit-cell along the ^-direction. 
The unit cell doubles because 5 = 1/2 models behave as if 
there is a topological flux piercing each unit cell, so (as for 
a magnetic field) the net topological flux through the cylin¬ 
der changes along the direction x. |28j Strictly speaking it is 
more precise to view the MES as a torsor for the fusion group, 
but this subtlety does not affect the measurements proposed in 
this work, so for notational simplicity we will label the MES 
by anyon types. 

III. SPACE-GROUP QUANTUM NUMBERS: ROBUST SET 
INVARIANTS 

In this section we consider the space-group quantum num¬ 
bers of a finite length cylinder: even though the structure of 
the two edges is non-universal, we argue that the space-group 
quantum numbers are. In particular, consider creating a pair 
of well-separated excitations from the vacuum which are re¬ 
lated by a mirror plane, and separating them out to the edges 
of the cylinder. If the reflection quantum number flips sign 
after this process, the excitations must be anyons which are 
connected by an invisible string which is odd under reflec¬ 
tion. Since the pair of anyons together transforms as R = — 1, 
it as if R = \J -1 acting on each anyon individually, so we 
say R is ‘fractionalized.’ 

While global quantum numbers are only well defined for fi¬ 
nite systems, at a later point we will show that they leave their 
imprint on the bulk entanglement spectrum in a way which 
can be measured on an infinitely long cylinder as well. 

A. Reflection quantum numbers 

Consider a large but finite cylinder with a reflection symme¬ 
try R x that exchanges its two edges. We argue that the global 
R x quantum number of any symmetric state | A, a) with no ex¬ 
citations in the interior of the cylinder depends only on a) the 
SET order of the bulk phase of matter b) the details of the ge¬ 
ometry, such as its dimensions, which we denote by ‘A’, and 
c) the topological sector a of the cylinder. Throughout this 
paper we will denote these global quantum numbers by Q, 

R x |A,a) = Qr x (A,o) |A,a). (1) 

To elaborate on c), note that an infinitely long cylinder has the 
same topological ground state degeneracy as the torus. While 
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the edges may reduce the ground state degeneracy, we only 
require that there are no excitations in the bulk of the cylinder, 
so are left with the same bulk degeneracy as the torus. We 
label these topological sectors of the finite cylinder by a, and 
as discussed earlier we assume a indexes a special ‘minimally 
entangled’ basis which has definite topological flux a thread¬ 
ing the cylinder. This work focuses mainly on Z 2 topological 
order because the states |A, a) will break R x if a is not it’s 
own anti-particle. 

The global quantum number Qr x is insensitive to any de¬ 
tails of the edge state or bulk Hamiltonian. To show insen¬ 
sitivity to the edge, note that perturbing the edges amounts 
to acting with unitaries Ul, Ur localized at the edges of the 
cylinder. When Ul and Ur are. spatially well separated, using 
R x symmetry we can require that Ur = R x UlR x 1 . Since the 
perturbation UlRxUlR x 1 commutes with R x , the quantum 
number is unchanged. Q is insensitive to the bulk Hamilto¬ 
nian because Q = ±1, so is quantized and can only change 
during a bulk phase transition. 

However, Qr x (A, a) can depend on the topological sector 
a, since changing the topological sector from a —>■ 6 • a re¬ 
quires separating an anyon pair b/b out to the edges using R x , 
which is a string-like operation. We will find that the depen¬ 
dence on A cancels if look at the relative quantum number 
between topological sectors, [57;] 


^.(b) _ Qrx (A, ba) 
Qfl.(A,a) ' 


( 2 ) 


These ratios have a particularly simple relationship to the SET 
order: if = — 1 , it implies that a pair of anyons related by 
R x each carry half the R x = — 1 quantum number, which we 
consider to be ‘fractional.’ In contrast, our earlier argument 
implies that a pair of truly local excitations must always have 

Qr* = 1 - 

Now suppose there is an additional reflection symmetry R y 
which does not exchange the two edges. In the absence of the 
edge-exchanging R x symmetry, the R y quantum numbers are 
not robust, because nothing then prevents the perturbation Ul 
from being odd under R y while Ur is even. But if both R x 
and R y are present, we can instead measure the combination 
I = R x R y \ since I exchanges the edges, according to our 
earlier reasoning Qj is also a protected invariant. 

Finally, if lattice doesn’t have a C 4 symmetry (as for the 
Kagome model) there are distinct ways to compactify the ge¬ 
ometry into a cylinder: for one choice, R x exchanges the 
edges, while for the second, R y does. We can then measure 
quantum numbers Qr x in the first cylinder, Qr v under the 
second cylinder, and Qj in either. This give three independent 
quantum numbers for each anyon type, Q^r^r which we 
will find almost fully characterizes Z 2 SETs (at least within 
the PSG framework). The remaining information relates to 
the commutation relations of time-reversal T and the reflec¬ 
tions R, which will lead to protected edge degeneracies we 
discuss in Sec lV Cl 


B. Translation quantum numbers 

In a magnetic field the translations T x , T y form a ‘magnetic 
algebra’ T x T y T x x T~ x = e 1 ® which is a projective repre¬ 
sentation of the translation group. Even in the absence of a 
physical magnetic field, in a topologically ordered phase the 
anyons may experience an effective magnetic field. The mag¬ 
netic field experienced by anyon a is encoded in the projective 
relation {T x T y T~ l T^)V> = r$. 

For an Abelian theory the projective relations must obey 
the fusion rule ri xy r] xy = r]£y\ since r/xj is the Berry phase 
acquired when a circles a unit cell. For a Z 2 spin-liquid in 
an S' = 1/2 model, we always have the relation r/xy = — 1. 
This can be argued in the language of Z 2 gauge theory, where 
all S = 1/2 objects are the source of Z 2 electric flux (for 
example the spinons b, /, which map on to the electrically 
charged e, / particles in the gauge theory). This includes the 
microscopic S = 1/2 in each unit cell, which implies that the 
system behaves as if there is electric flux piercing each unit 
cell. Consequently under e-m duality, the flux m (the vison) 
experiences a background flux of 7r per unit cell, so r ] xy = 

-1. Since r^y = l,rji v J = -1, and rjxyvi'y = viy\ there is 
a single sign left undetermined, which is the most basic SET 
distinction between Z 2 spin-liquids. 

To probe r/ xy , consider a cylinder of length L x and circum¬ 
ference L y in topological sector a, and measure the momen¬ 
tum quantum number Qt v (A, a) = T y \ A, a). Note that when 
L y is odd the MES double the unit cell in the x direction, 
so we must restrict to L y even. T y symmetry alone does not 
protect Qt v , since the edge excitations can carry an arbitrary 
momentum, but the combination of T x and T y allows us to 
define a robust ‘momentum per unit length.’ Recall that in 
the Landau gauge, the momentum T y of a particle in a mag¬ 
netic field is proportional to its position x. Since | A, a) has 
an a particle localized near the edge, as we grow L x its mo¬ 
mentum grows linearly with x, ie, it is a momentum per unit 
length of cylinder. To define the momentum per unit length 
operationally, we need to grow the length of the cylinder L x 
while keeping the topological flux and edge state the same, 
as otherwise Qt v could contain a spurious contribution com¬ 
ing from the changing edge state. Concretely, we require the 
reduced density matrix for the edge be kept constant as L x 
grows. The momentum per unit length r/xy is then 

Qt v (A, a) = q{dA, a) . (3) 

q(d A, a) depends on the edge, while the bulk contribution re¬ 
veals the SET invariant r]xy . 

The momentum per unit length is trivial to measure in any 
tensor network ansatz. For MPS, it is the invariant r/xy = 
VT y defined in Eq. (l48l) i when viewing g = T y as an ‘onsite’ 
symmetry in the ID representation of the cylinder. tt v is a 
byproduct of the algorithm used to calculate the momentum- 
resolved entanglement spectrum of infinite-DMRG studies, so 
presumably has already been computed in existing studies. 
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IV. SYMMETRY ENRICHED ORDER: DETECTING THE 
PROJECTIVE SYMMETRY GROUP 

The preceding discussion is independent of any classifica¬ 
tion of space-group SETs, since we have argued on general 
grounds that these quantum numbers are robust SET invari¬ 
ants. Nevertheless, we would like to identify these invari¬ 
ants within a general classification scheme. Currently there 
is not a complete classification of space-group SETs. How¬ 
ever, the parton construction provides a rich zoo of Z 2 spin- 
liquids whose symmetry properties can be analyzed using 
Wen’s ‘projective symmetry group’ (PSG). 0 The PSG pro¬ 
vides at least a partial classification of space-group SETs. In 
this section, we show how to compute the quantum numbers 
Qu (A, a) within the parton construction, thereby identifying 
them with invariants of the PSG. 


A. The Parton Construction and the Projective Symmetry 
Group 

The resonating valence bond (RVB) picture proposed by 
Phil Anderson provided the first intuition for a spin liquid 
ground state as a quantum superposition of different dimer 
configurations covering a lattice of spin-1/2 particles. Each 
dimer (denoted by • — •) is a singlet pair formed by two spin- 
1/2 particles: 

!•-•} = | t)lU)2-U)l|t)2. 

The state is a ‘liquid’ because the quantum superposition of 
dimer patterns restores the translational symmetry. One type 
of elementary excitation in these systems is created by break¬ 
ing a dimer into a pair of particles carrying spin-1/2 each, 
which were coined spinons. 

The parton construction is a systematic formalism for writ¬ 
ing down ansatz RVB wavefunctions in which each spinon is 
realized either as a fermionic parton f a or bosonic parton b a , 
where 0 =}, }. In the fermionic description, each dimer is re¬ 
alized as an s-wave Cooper pair of partons; breaking a Cooper 
pair generates a pair of spinons. The microscopic spins S r are 
related to the partons through the bilinears 

‘“’r = 9 y ] fr,a a ot,pfr,P = ^ ^ ' ^]-,a (J a,/3^r,/3-(4) 

a,0=t,l 

where o are the three Pauli matrices. In order for this mapping 
to generate a sensible S = 1/2 wavefunction, the partons can¬ 
not be free particles: they obey the “single-occupancy” con¬ 
straint of one parton per lattice site: 

X fr,afr,a= X t = 1, V lattice site r. (5) 

(7=1,1 0-=t,l 

This constraint can be implemented by a gauge field which 
couples to the partons. 

In practice, we use the parton construction to create ansatz 
wave-functions. If |MF) is a ‘mean-field ansatz’ state for the 
partons which need not obey the single occupancy constraint 
(for example, a BCS superconductor of fermionic partons f a ). 


we enforce the constraint via Gutzwiller projection to obtain 
an S' = 1/2 wavefunction: 

(ti +2 = 

(Ol/r^l/r^v/r.^-'-lMF}, (6) 

and similarly in the bosonic construction. Note that in the 
fermionic case, we must choose and fix an ordering of sites 
r i, r 2 , ■ ■ • in order to maintain the correct relative sign be¬ 
tween different spin configurations. For the purposes of cal¬ 
culation |MF) is usually taken to be a free wavefunction, such 
as a mean-field BCS superconductor or pair-superfluid for the 
fermionic / bosonic constructions considered here. Gutzwiller 
projecting the creation of a single parton f'/b' results in a 
highly non-trivial excitation: an S' = 1/2 anyonic excitation, 
the spinon. 

A crucial question in the parton construction is how the 
symmetries {U £ SG} of the 5=1/2 wavefunction (such 
as global SO(3) spin rotations and space-group symmetries) 
are realized in the partons and their mean-field ansatz |MF). 
The simplest possibility is that the partons form a linear rep¬ 
resentation of the symmetry group SG. But is this the only 
option? The answer is no, because according to Eq. iQ} apply¬ 
ing a gauge transformation b rrr —y e'^ r b r a leaves the physi¬ 
cal spin operators unchanged, and accordingly the Gutzwiller 
projection in Eq. (| 6 ]) is a many-to-one mapping insensitive to 
U( 1) gauge transformations. Consequently under a series of 
symmetry operations {L/} £ SG which yield the identity op¬ 
eration e 

U 1 U 2 ---U n = e ( 7 ) 

a single-parton operator b r (T (or /) may acquire a nontriv¬ 
ial phase factor e 1 / 1 instead of remaining invariant. In 
this case the partons transform projectively, rather than lin¬ 
early, under the symmetry group SG. The symmetry oper¬ 
ations {U £ SG} are accompanied by certain gauge trans¬ 
formations {Gu\U £ SG} on partons, forming a “projec¬ 
tive symmetry group” PSG = {GuU\U £ SG} which is a 
central extension of the symmetry group SG. [T|] The cen¬ 
ter of such an extension is called the “invariant gauge group” 
IGG = PSG/SG. The IGG are those gauge transformations 
which leave the mean-field ansatz |MF) invariant. 

For the Z 2 spin liquids which are the focus of this work, the 
invariant gauge group is IGG - Z 2 . This means under a series 
of symmetry operations {(/,} in ([7}. each parton acquires a 
Z 2 phase factor of ±1. This is consistent with the mean-field 
ansatz being a singlet BCS superconductor of spinons. For 
example, for symmetry a symmetry UU = e we have the 
relation 

(GuUffr^iGuU)- 2 = V f u f ri , ai , rb = ± 1 . ( 8 ) 

The primary goal of this work is to show how to detect these 
±1 phases associated with symmetry group SG. 

B. Minimally entangled states within the parton construction 

We first must review how to generate the MES within the 
parton construction. 
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There is a two-fold degeneracy associated with threading 
a vison v for both finite and infinite cylinders. In the parton 
ansatz this arises because the boundary conditions of a cylin¬ 
der of circumference L y can be either periodic (P) or anti- 
periodic (AP). At the level of the parton Hamiltonian, this is 
accomplished by assigning an additional sign of —1 to all ma¬ 
trix elements which cross a line at some fixed y = yo', the 
choice of yo is a gauge choice. When acting with T y or R y , 
the location of the twist yo must be restored by an additional 
contribution to the gauge transformations Gt v /r v , leading to 
new PSG relations: 

{G T T y ) L * = (-l) bc 
(9) 

(GT y Ty) Ly / 2 (GR v Ry)(G Ty Ty)~ Lv ^ 2 (G Ry Ry)~ 1 = (~1)^ 

( 10 ) 

where be = 0/1 for PjAP. Hence in all that follows the 
PSG implicitly depends on the boundary condition (be) of the 
mean-field ansatz |MF) 

The two-fold degeneracy associated with threading a 
spinon is a bit more subtle, as in the finite case the degeneracy 
is split by the edges. If we make a bipartition of the cylin¬ 
der at some xq, the parton parity (—1 ) N Wf in the left half of 
the mean-field ansatz |MF) fluctuates across the cut (note the 
parton number itself is not conserved). After Gutzwiller pro¬ 
jection, Eq. <(6j). the parton parity to the left is fixed by the 
number of sites to the left. But in the infinite case there is an 
ambiguity, since the number of sites is infinite. This means 
that when Gutzwiller projecting we can freely choose either 
the sector with even (E) or odd (O) parton parity to the left of 
the cut at Xq, which generates an additional 2-fold degeneracy 
on the infinite cylinder. 

The choice of P / AP boundary conditions combined with E 
/ O parton parity generates the 4-fold degeneracy of the infi¬ 
nite cylinder. These sectors are identified with the anyon types 
in a manner that depends on the parton construction: 

1, (P, E), (AP, E), (P, O), (AP, O) (bosonic) (11) 

1, v,b,f<-> (AP, E), (P, E), (AP, O), (P, O) (fermionic) 

( 12 ) 

Note the role of P/AP is flipped between the two construc¬ 
tions. 

For an even circumference cylinder the E / O parity is the 
same for all cuts Xq, since an even number of sites intervene 
between cuts. But for an odd circumference cylinder, the E / 
O assignment alternates with xo ■ This alternation doubles the 
physical unit cell. 

C. Computation of global quantum numbers from the PSG 

Since the PSG determines how the partons transform under 
symmetry operations, we can compute the (crystal) symme¬ 
try quantum numbers of any projected wavefunction (Eq. ([6])) 
constructed from a parton mean-field ansatz. If U is a space- 
group operation which permutes the sites according to i —> 


U(i), the wavefunction transforms as 

= <W- 1(i) }|*> = <0|nw ito IMF), 

(13) 

and similarly for the bosonic ansatz. We can split the resulting 
quantum number into two contributions. First, there is a part 
Qu (A, be) coming from the Gutzwiller projection, 

< o in^^- 1(i )=^( A > bc )< o in/^^ 0 4 ) 

i i 

which depends only on the geometry A and the PSG (which 
is modified by be); we will show how to compute this in the 
subsequent section. The second contribution comes from the 
quantum number of |MF). Inserting Eq. (1 1 4b into Eq. ( 1131 ) we 
find that the total quantum number factorizes as 

Gull |MF) = Qtr(MF) |MF) (15) 

U |d>) = Qj/(A,bc)Qt/(MF) |\P) (16) 

D. Ratios of edge-exchanging quantum numbers in different 
topological sectors 

In the first scenario, we suppose we have access to all sev¬ 
eral topological sectors on the same geometry. We find that 
for edge-exchanging symmetries, the ratio between the quan¬ 
tum number before and after threading anyonic flux a reveals 
the PSG of anyon a. 

1. Spinon insertion 

We fix the geometry A and compute the relative reflection 
quantum number Qr between two states which differ by the 
insertion of a pair of spinons at the edges. To generate the 
appropriate pair of mean-field ansatz |MF) and |c • MF) which 
differ by spinon flux c = b/ f (depending on the construction), 
let 4 create an arbitrary bosonic / fermionic parton near the 
left edge. To ensure that c • MF is symmetric under R we 
must create a corresponding spinon on the right using 4 = 
(GrR)cI(G r R)-\ so 

|c • MF) = 4 (GrR^GrR)- 1 |MF). (17) 

Qu{ A, be) is unchanged, and it is straightforward to verify 

GrR\ b/f ■ MF) = (-1) f (G r R) 2 Qr(MF) \b/f ■ MF) 

(18) 

where the sign (—1) F occurs forfermionic partons as we must 
exchange the creation operators and we use (GrR) 2 to denote 
the parton PSG associated with R 2 = e: 

(GrR) 2 Ci(G r R)~ 2 = rift Ci, V*, 

(GrR) 2 = rf k . 


( 19 ) 
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Consequently for any geometry the change in the quantum 
number Qr after inserting a bosonic / fermionic spinon is 


n (b/f) = Q fl (A,b//-MF) 

Qr(A,MF) 




( 20 ) 


Hence the spinon PSG can be recovered by measuring the rel¬ 
ative quantum number between different topological sectors. 

The relative quantum numbers can be computed for any 
space-group symmetry which exchanges the edges; by us¬ 
ing different cylinder compactifications we can measure the 
R x ,R y , and I = R x R y quantum numbers. There may be 
distinct 7r-rotations depending on whether the rotation is site 
or bond / plaquette centered; on a square lattice, for example, 
/' = T x R x T y Ry, reveals an independent PSG relation. In 
Figs. m and UH we tabulate the relative quantum numbers for 
the square and Kagome lattices. 


2. Vison insertion 

We again fix A, and compute the relative reflection quantum 
number Qr between two states which differ by the insertion 
of a vison. As discussed, threading a vison switches between 
P/AP boundary conditions, so Qu{ A, be) may change due to 
the PSG relations of Eq. <(9]). The PSG relations associated 
with U 2 = e are only modified if U takes an odd number 
of sites across the twist boundary condition modified by the 
vison. For both the square and Kagome lattice, for the ge¬ 
ometries in which U exchanges the edges of the cylinder we 
have 


Q flH (A,AP) = Q fly (A,P) 

(21) 

Qr*(A,AP) = Qr x (A,P) 

(22) 

Q /6/p (A,AP) = Q 7b/p (A,P) 

(23) 

Qj, (A, AP) = -Qj, (A, P) 

(24) 


where I b / p is a is bond or plaquette centered 7r-rotation, while 
I s it site-centered. 

Next we argue that Q[/(MF) is unchanged when threading 
a vison. To start, suppose |MF) is in the same phase as the 
ground state of a BCS superconductor / pair super-fluid. In 
the subsequent section, we show that <5 j/(MF) = 1 regardless 
of boundary condition, so is unchanged by vison insertion. 
Now suppose we modify the ground state |MF) with arbitrary 
[/-symmetric edge perturbations of fixed parton parity. The 
resulting Q[/(MF) can depend only on the parton parity of the 
edge perturbation, but not on the boundary condition, because 
the vison-modified PSGs of Eq. (0 will necessarily act on U- 
related partons and the signs will cancel. This shows Qij(MF) 
is unchanged by vison insertion, and we conclude that 


Qr! = 1, 


Q 


( v ) _ 


= 1 , 


QrI = 


q 


=-i. 
(25) 


A possible loophole in our argument is that we assumed that 
for the ground state, |MF) could be taken to be in the same 
phase as a BCS / pair superfluid, but our result agrees with 
earlier discussions. ||5ll 


a 



b 



f 

V 

= 

Q (v 

Q (b) 

QrI 



VR X 



~Vr x 

1 

Vr x = 

-Vr 





VRy 



-VLy 

1 

VRy = 

~Vr 



VR x 

VR y 

VR X ,Ry 

-4, 

VLy 

^RtcRy 

1 

Vr x ,r v 

= - 

7 lR x ,R y 


vlyVR x 

VR y 

Rx iRy 

-vlyVR, 

VLy 

^ R X ,Ry 

-1 

Vxy 

~Vxy 



TABLE I. Relative quantum numbers between topological sec¬ 
tors a of a square lattice spin-liquid. I p / b is a plaquette or bond 
centered 7r-rotation, while I s is a site-centered 7r-rotation. 


a 

b 

/ 

V 

q(/) = qWqW 

eg 

(_1)P2+P3 

-V* 

1 

(_ 1)P 2+P3 = _ Va 

eg 

(_1)P2 

~Vcr Vo-Cq 

1 

(-1) P3 = Vo Ce 

Qa 

(_1) P1 + P3 

~VC e 

1 

(_l)Pi = -•q Ca r) aCe 

Q V 

(-!) p 3 

-V12 VC 6 

-1 

(_l) p i = - Vl2 


TABLE II. Relative quantum numbers Q between topological sec¬ 
tors a in a Kagome lattice spin-liquid. Ih/ S is a hexagon / site cen¬ 
tered 7r-rotation. The bosonic PSGs are expressed through the invari¬ 
ants (pi,p 2 ,P 3 ) of Ref. & while the fermionic PSGs are expressed 
through the invariants r of Ref. m . The fermionic invariants satisfy 
Vr2VC 6 VaC 6 = 1 tautologically. 


3. Fusion and unification: Q^R = Q60q(&) 

Our derivation of the vison quantum numbers shows that 

( f) 

the relative reflection quantum numbers obey fusion, = 
Qu^Qu *’ t> ecause vison insertion changes the quantum num¬ 
ber by regardless of the parton parity at the edge. As 
tabulated in Table U and Table HT1 we can use this fusion rule 
to equate the bosonic and fermionic PSGs, unifying the two 
approaches. 


E. The quantum numbers of the ground state 


A given geometry A genetically has a lowest energy state 
which is SO(3) symmetric on both edges, and two-fold degen¬ 
eracy for vison insertion. We can compute the quantum num¬ 
bers for such an SO(3) symmetric state by assuming |MF) is 
in the same universality class as a BCS superconductor / pair 
superfluid in the fermionic / bosonic constructions. |MF) is 
then invariant under any symmetry: Qc/fMF) = 1. This is 
because in both constructions we have 


|MF) = exp [Eij ff (» - i)4Al\ i°}’ 
ff(i-j) = (—l) F 3(j -i) 


(26) 


where c = f/b and (— 1) F = —1/1 for Abrikosov-fermions 
and Schwinger-bosons respectively. Because |MF) always 
contains the parton Fock vacuum |0) as a component in the 
Taylor expansion of the exponential, and |0) is neutral under 
any symmetry operation, |MF) must also be neutral. There¬ 
fore the quantum number of the ground state depends only on 
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the geometry A and the PSG: 

Qt/(A,bc,MF) = Qu(A,bc) (27) 

where be will depend on the two-fold degeneracy associated 
vison insertion and Qu(A 1 be) is defined in Eq.[l4] This result 
is particularly useful numerically, since by modifying A we 
can probe the PSG using only a single topological sector. 

In the following we specifically illustrate how to obtain the 
eigenvalues Qi;(A. be) of two crystal symmetry operators, in¬ 
version I and mirror reflection R, for a projected wavefunc- 
tion iona finite-size lattice from parton PSGs. 

1. Eigenvalue of plaquette-centered inversion I p 

First we consider an inversion symmetry I p whose inver¬ 
sion center lies on a plaquette. For a finite-size lattice with 
plaquette-centered inversion I p , the number of lattice sites N s 
must be even. All lattice sites must be exchanged in pairs 
under inversion operation, since no lattice site remains invari¬ 
ant under I p operation. More specifically, spinon f ri , ai must 
appear altogether with its inversion counterpart ff ra = 

Gi p Ip / r , ,ai(Gr p Ip)- 1 in the many-spinon operator J], /r^- 
Note that in projected wavefunction ((gJ there is always a par¬ 
ticular ordering for the real-space positions {iy} of the many- 
spinon operator fl,; ■ Here we simply choose a ordering 
in which a pair of spinons related by inversion show up to¬ 
gether i.e. 

n, fw = n' fw • G Ip i p f r ,. a , {G Ip i p )-\ ( 28 ) 

where ' denotes the product over half lattice sites that are un¬ 
related by inversion. Clearly under inversion operation I p the 
above many-spinon operator transform as 

GiJpUJ^MGiJp)- 1 = 
U'GiJpfr.MGiJp)- 1 ■ (.GiJpYfr^AGiJp)- 2 . 

By definition of PSGs we have 

( Gi p Ip) fn^iiGipip) = Vi p fri,cTii Vi p = ± 1-(29) 

since I p = e yields the identity operation. We used rj- to 
denote the PSGs for fermionic spinons (Abrikosov fermions) 
and rj b for bosonic spinons (Schwinger bosons). Notice 
that for Abrikosov-fermion representation, exchange of two 
spinons f Ti and /) r gives rise to an extra —1 sign due to 
Fermi statistics. As a result we obtain 

GipIpUfr^iGipIp)- 1 = (-V f Ip ) Na/2 Uf^ ( 3 °) 

i i 

Hence the eigenvalue of plaquette-centered inversion I p for 
projected wavefunction ([6} on a Ay-site lattice (N s =even) is 

Qi P (N s ) = (-v {) N ‘ /2 (31) 

for Abrikosov-fermion representation. On the other hand, in 
the Schwinger-boson representation, exchange of two spinons 
won’t yield a —1 sign and we have 

Qi P (N s ) = ( V b Ip ) Ns/2 (32) 


2 . Eigenvalue of site-centered inversion I s 

Now let’s take one more step to consider an inversion sym¬ 
metry I s whose inversion center lies on one or more lattice 
sites. Fet’s assume inversion centers contain Nj sites and N s 
is the total number of lattice sites. For those (N s — Ni) sites 
other than the inversion centers, their contribution to the I s 
eigenvalue follows exactly the same form as ( |3T1 ) and ( 132k 
except that we need to replace N s by N s — Nj. 

What about the contribution from the Nj inversion centers? 

First of all, if there is an odd number of inversion centers 
(Nj =odd), the inversion eigenvalue is not a gauge invari¬ 
ant quantity since the symmetry operations on a single spinon 
can always be followed by an arbitrary gauge transformation. 

In the case when Nj =even, if the inversion centers are not 
related by any other symmetry, again they can acquire extra 
gauge transformations independently under symmetry opera¬ 
tions, and again the inversion eigenvalue is not a topological 
invariant. 

If an even number of inversion centers are related by sym¬ 
metry, on the other hand, one can compute their contribution 
to I s eigenvalue from parton PSGs in a universal manner. 
Without loss of generality, let’s consider a pair of inversion 
centers related by certain crystal symmetry P (e.g. it could 
be a mirror reflection or a translation on a finite cylinder), this 
spinon pair operator transforms under inversion I s as 

Gljs [f T ,a ■ GpPfrAGpP)- 1 ] (Gish )- 1 
= GjsIsfrAGish )- 1 ■ GjJApPfrAGpPrYGjsis )- 1 
= vi^GishfrAGish )- 1 ■ GpPGjsisfrAGisQ-YGpP )- 1 

= vi, P vi ■ GpPfrAGpP)- 1 ] ■ (33) 


GisisGpPfrAGpPrpGjsis)- 1 = (34) 

V f Is , P GpPGjshfrAGishrHGpP)- 1 , 

(GlsIsYfrAGlsh)- 2 = Vis /r,T- (35) 

Consequently, for a Aj s -site lattice with Ay inversion cen¬ 
ters ( N s ,Ni =even) which are pairwise related by crystal 
symmetry P, the inversion eigenvalue of projected wavefunc¬ 
tion ([6]) is given in terms of spinon PSGs by 

Qis(N s ,Nj) = (-v I i a ) ( ‘ N '~ Nx)/2 ( , n f i a v I i t „p) Nl/2 

= (-V f is ) Ns/2 (-v{,p) Nl/2 - (36) 


Qi e (N s ,Nj) = (v b j s ) N ‘ /2 (v b is,p) Nl/2 - < 3? ) 
for Schwinger bosons. 

A crucial point is that the PSG rji e ,P can depend on the 
boundary condition of the cylinder (for example, if P = 

Ty' 2 ). 


where we defined spinon PSGs 


for Abrikosov fermions and 
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3. Eigenvalue of mirror reflection operator R 

The eigenvalues of mirror reflection operator R can be com¬ 
puted completely in parallel to the case of inversion symmetry 
as discussed previously. Again let’s assume Nr lattice sites 
lie on the mirror reflection axis on a A\ s -site lattice. As argued 
earlier, only when Nr is even and these Nr sites are related to 
each other by other crystal symmetries, will the R eigenvalue 
be a topological invariant that is fully determined by parton 
PSGs. Let’s assume these Nr sites are exchanged in pairs by 
crystal symmetry P. Similar to the case of inversion symme¬ 
try I s we can compute the reflection eigenvalue of projected 


wavefunction ([6]) as 

Qr(N S: Nr ) = {-rf R ) N - /2 {-rf R<P ) N « /2 - (38) 

for Abrikosov fermions and 

Qr(N s ,Nr) = (v b R ) Nj2 (v b R,p) NR/2 - (39) 

for Schwinger bosons. The parton PSGs are defined as 

G r RG p P f r ,o{GpP)~ 1 (G R R)- 1 = (40) 

V ^ P GpPGRRf r , a {G R R)-\GpP)-\ 

0 G R R) 2 frAG R R)" 2 = v f R f P *. (41) 


for Abrikosov fermions and similarly for Schwinger bosons. 


4. Unifying bosonic and fermionic PSGs 

If an Abrikosov-fermion state and a Schwinger-boson state 
describes the same Z 2 spin liquid state, their symmetry quan¬ 
tum numbers on any finite lattice must be the same for ar¬ 
bitrary crystal symmetries. Therefore from the eigenvalues 
of inversion and reflection symmetries summarized previ¬ 
ously, we can achieve a unification of Abrikosov-fermion and 
Schwinger-boson representation: i.e. their PSGs must satisfy 


the following correspondence: 

~ V f i p = V b i p ', (42) 

-V f i e =Vif, (43) 

-vi,P = V b i s , P , (44) 

V crystal symmetry P satisfying PI S = I S P; 

-Vr = r) b R ; (45) 

- ^r^p = t Ir,p^ ( 46 ) 

V crystal symmetry P satisfying PR = RP. 


These relations are in agreement with our conclusions based 
on the relative quantum numbers. In the next section we’ll 
establish the correspondence between Schwinger-boson and 
Abrikosov-fermion representations for those PSGs concern¬ 
ing time reversal symmetry. This is achieved by relating the 
2D parton PSGs to ID SPT invariants by considering pro¬ 
jected wavefunctions on a thin but long cylinder. 


V. DIMENSIONAL REDUCTION AND ENTANGLEMENT 
SIGNATURES 

A 2D model defined on a cylinder can be viewed as a ID 
system by grouping one ring of the cylinder into a single 
super-site. This point of view is useful because the interplay 
of symmetry, topology and entanglement has been completely 
understood in ID through the recent classification of ID SPT 
phases. EMI In this section we explain how 2D SET or¬ 
der manifests itself as ID SPT order under this dimensional 
reduction. In particular, we find the Z 2 PSG relations have 
a one-to-one correspondence with the U(l) projective repre¬ 
sentations that classify ID SPT phases. While generally U(l) 
projective representations are a coarse-grained version of Z 2 
projective representations, space-group symmetries actually 
have an anti -unitary character under the dimensional reduc¬ 
tion, and for this special case the correspondence becomes 
one-to-one. 


A. A review of ID SPT phases. 

While the classification of ID SPTs can be discussed in 
terms of Schmidt decomposition, the most compact treatment 
uses the formalism of matrix product states. We refer to pre¬ 
vious works for a more detailed review. EMI 


1. Matrix product states 

Let | j n ) span the local Hilbert spaces of a spin chain with 
sites at n. A MPS I’L) is characterized by a sequence of rank- 
3 tensors {T^" ail+1 } and rank-1 vectors {s a „} through the 
ansatz 

Xn 

E ( 47 > 

ot n — 1 n 

The indices a n which are summed over are called the auxil¬ 
iary indices, with dimension \n- We have assumed the MPS 
is in the ‘canonical form,’ which means that each s a „ is the 
set of Schmidt weights for a bipartiton of the system between 
sites n — 1 ,n. 

The MPS ansatz includes both finite and infinite spin 
chains. In the finite case with L sites, xi = Xl+ 1 = 1. In the 
infinite case with a unit cell of length L, we can always choose 
the tensors to share this unit cell: T 3 u p p LC « n+1+i = an+1 
and likewise for s and X- 


2. Onsite symmetries 

If a spin chain is invariant under an onsite symmetry g £ G 
(e.g. a spin rotation) , it is natural to ask how the symmetry 
is encoded in the tensors T. The representation of the onsite 
symmetry decomposes into its action on each site, g = <E> n g n - 
For notational simplicity, we will drop the site index n. An 
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MPS is symmetric under g if and only if the V and s transform 
as m 


E 9A =r g J2 UsA'Ap ( 48 > 

k a',/3' 

/3 — S a Ug; a p (49) 

where the U g are unitary matrices and r g £ U(l). 

The phases r g form a U(l) representation of the group, and 
encode the g-charge per unit length. But for the unitaries Uh 
there is another possibility: U g may be a projective represen¬ 
tation, meaning that the requirements of a group representa¬ 
tion are satisfied only up to U(l) phases. This subtlety arises 
because the phase of U g isn’t fixed by the transformation law 
of Eq. @ 8 ] Arbitrarily fixing the phase of each U g , the U(l) 
phases to are encoded in the relations 


There is also an interplay between the other symmetries and 
T. If G contains H x Zj C G as a subgroup, where H is 
onsite, for h £ H there is a projective relation 


UtKUt 1 


u(r,h) TT 

u(h,T) h 


(56) 


For our purposes it will be sufficient to understand A = Z 2 ; 
since aT is itself a unitary Z 2 symmetry we can form the 
gauge-invariant relation 


UhtUht = 7hr = ±1, (when h 2 = 1), (57) 

an additional Z 2 invariant. Alternatively, we have the relation 

{UlWrKUrA) = IhTlr (58) 


U g U h = u{g, h)U gh , u>(g, h ) £ U(l). (50) 

The phases ui are called cocycles, or the factor set. If we con¬ 
sider the phase ambiguity U g 6(g)U g , 9(g) £U(1) to be a 
‘gauge transformation,’ we see that ui is not gauge invariant. 
The classification of gauge-inequivalent ui is given by the 2nd 
group cohomology [w] £ 7f 2 [G, 17(1)], resulting in a classifi¬ 
cation of ID phases symmetric under G. Note that for sym¬ 
metry group G = Zjv, there are no projective representations, 
and hence no ID SPT phases. 

An important physical signature of an SPT phase under on¬ 
site G is degenerate edge states. It can be shown there are 
edge states which transform under G with the same projective 
representation [w] as the U g \ if [w] is non-trivial, the projective 
representation must be multi-dimensional, implying a degen¬ 
eracy. 


3. Time-reversal symmetry 

The transformation laws are modified for time-reversal 
T = ® n urK because of the complex conjugation K. Similar 
to before, the MPS can be taken to transform as 

E U T;A = E UwK'A'P ( 51 > 

k <x' ,/3' 

Ur-,apsp = s a Ur- a p (52) 

(note we can remove rj- by redefining the U(l) phase of the 
state). But for an anti-unitary symmetry like time-reversal T, 
the projective relations are modified to 


4. Reflection symmetry 

Finally, consider a reflection R which spatially inverts the 
ID chain, R \j n ) = UR-j nt k_„ \k- n ). The unitary matrix ur 
encodes any internal rotation in the definition of the reflection. 
The transformation law is 0 

E E = r R E (^9) 

k ex' ,/3' 

UR, a pSp = SaUR-ap (60) 

where vr = ±1 provides the first Z 2 invariant, the ‘parity per 
unit length.’ Somewhat surprisingly, combining this transfor¬ 
mation law with those of an onsite h we find the projective 
relations are the same as those of a anh-unitary symmetry: 

U R Uj = uj(R, h)U Rh , (61) 

U h U R = u(h, R)U hR u £ U(l). (62) 

The origin of this similarity is that transposition T and com¬ 
plex conjugation behave analogously when acting on the uni¬ 
tary Uh- This point is important, as it implies inversion has 
the same Z 2 invariants as time-reversal T: 

UrUr = 7 a = ±1 (63) 

VhrUIr = ihr = ±1 (when h 2 = 1) (64) 
{U 2 h )(U R U* h U^U^) = lhRlR = ±1 (65) 

B. Identification of ID SPT order and the space-group PSGs 


U T Uh = w(T, h)U T h, (53) 

U h U T = u(h, T)U hr to £ U(l) (54) 

In contrast to an onsite G = Z 2 , for the anti-unitary time- 
reversal (G = Z 2 ') we can form the gauge-invariant relation 

7r = U T U* r = ±1. (55) 


Earlier we argued that space-group quantum numbers are 
topological invariants in the presence of reflection symme¬ 
tries, and calculated these quantum numbers using the PSG. 
We now show that under the dimensional reduction these 
quantum numbers can be calculated from the ID SPT invari¬ 
ants tr, 7 h, 'jhR- This clarifies the origin of their stability, 
since ID SPT phases are robust in the presence of symme¬ 
tries, and provides a dictionary between the 2D SET and ID 
SPT order. 


This 77 - gives a Z 2 SPT classification. 
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To equate the 2D SET order with the ID SPT invariants, 
we must compute the quantum numbers of a finite chain given 
the infinite chain ID SPT invariants. For notational simplicity, 
we will assume the dimensional reduction results in a one-site 
ID unit cell. The notation is more complex for the Kagome 
model, since under the dimensional reduction it is more con¬ 
venient to use a unit cell of two, so we delay the discussion 
of this case. Since the system is a 2D cylinder, in addition to 
a reflection R which exchanges the edges of the chain, there 
may be an orthogonal reflection R' which does not exchange 
the edges, but instead behaves as an ‘onsite’ symmetry under 
the dimensional reduction. We combine these two reflections 
to form a 180 degree inversion I = RR', which again ex¬ 
changes the edges. Consequently there are Z 2 ID SPT invari¬ 
ants = 7 RR'. 

For any R/I symmetric state on a finite chain with L sites 
(L-odd implies a site-centered inversion) it can be proven that, 
independent of any details of the edge, the quantum numbers 
of the chain are 

Qr{L) = 7 r ( r R ) L ( 66 ) 

Qi{L) = 7 / {ri) L (67) 

where 7 , r are the ID SPT invariants of the bulk phase. 

To relate the ID SPT invariants to the PSG, fix the trans¬ 
verse geometry Aj_ of the cylinder (such as the circum¬ 
ference), and find the four ground states of the infinite 
cylinder. Each topological sector a has ID SPT invariants 
7t/(A_i_, a), rr/(Aj_, a). 

Following our earlier discussion, we first compute the ratio 
of quantum numbers after threading a spinon b/f: 


Q (b/f) = Qu(A,b/f -a) = 7 t/(A_L, b/f ■ a) , rjj(A±, b/ f ■ a) 
u Qu(A,a) 7 c/(A_L,a) ny(Aj_,a) 

( 68 ) 


Since the result is independent of the geometry, we have the 
following: 


n (b/f) _ 7 t/(Aj_, b/f ■ a) 

MA±,a) 

ru(A±,b/f ■ a) = ru(A±,a) 


(69) 

(70) 


We find that the ratio of ID SPT invariants 7 jj between the 
topological sectors of an infinite cylinder reveals the PSG 
relation Furthermore, the parity per unit length is 

unchanged by threading a spinon. 

We then compute the ratio of quantum numbers after 
threading a vison v: 

0 (v) = Qu(A,v a ) = lujA^v ■ a) ru(A ± ,v ■ a) L 
u Qu(A, a) 7 £/(Aj _,a) ru(A ± ,a) 

(71) 


When U is not site-centered, our earlier discussion found that 
Qu' 1 = 1 , so by comparison we expect 

1 u(A_l,v ■ a) = 7 c/(Aj_, a) (72) 

ru(A±,v ■ a) = ru(Aj_,a), UG{R,I h / p }. (73) 


On the other hand, the site-centered inversion I s requires odd 
L, and from = — 1 we find 


r h ( A _l, v • a) 
r i b (A_l, a) 


(74) 


1. Translating between the Z 2 PSG and the ID SPT U(l) PSG 

The bosonic Z 2 PSG is encapsulated by the projective rela¬ 
tions 


X-Y = uj b {X,Y)(XY) (75) 

where u b € Z 2 . The ID SPT relations are similar, but 
u> € [7(1), which is in general a much coarser classification 
as there are more phase ambiguities. Yet we have shown that 
the ID SPT relations recover the 2D PSG. So how do the 2D 
PSG relations ‘descend’ to the ID SPT relations? The key 
point is that if the symmetry R exchanges the edges of the 
cylinder, we can transcribe the 2D PSG relations into ID SPT 
relations if we remember that reflection / inversion becomes 
anti-unitary in the ID SPT realization: 

R 2 = Vr=> U r U* r = lR (76) 

Due to the anti-unitary nature, the U(l) phase ambiguity does 
not affect the robustness of j R . In general, the 2D PSG 

relations descend to ID SPT relations if we treat edge¬ 
exchanging symmetries as anti-unitary. This explains the 
equality of the quantum numbers found in Eq. ( l69l ). 

) C. Identification of ID SPT order and the RTR~ 1 T ~ 1 PSGs 

In the fermionic parton construction there is an additional 
fermion PSG associated with the interplay of time-reversal T 
and a reflection R: 


r] f RT = R~ 1 T~ 1 RT (77) 

Viewing R as on-site, under the dimensional reduction we 
have a similar ID relation 


UrU t = ^ ^U T U* R (78) 

w(T, R) 

But lu(R,T) is not U(l) gauge invariant. Instead, we may 
consider the ID SPT invariant defined by Eq. (157k Comparing 
the 2D and ID SPT PSG relations, 

(R 2 ){R- 1 T~ 1 RT) = ri f R ri f RT 2D PSG (79) 

{U R ) 2 {U R 1 U T U* R Ufi 1 )= lRTl T ID SPT (80) 

we obtain the following 2D PSG to ID SPT reduction: 

VrVrt 1RT1T- (81) 

This identification can be verified by checking for the physical 
signature of t) R T) RT = —1. The fermionic spinon always has 
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a 2-fold degeneracy associated with it’s Kramer’s degeneracy 
77 j, = — 1 , but when 7 ^ 7 ^ = — 1 the spinon excitation is at 
least four-fold degenerate, as this is the minimal dimension of 
these projective relations. The ID SPT relation 7 t = —7 rt 
implies precisely this additional two-fold degeneracy. 

In summary, we can determine 7 /^ 7 ^, T by computing the 
change in the ID SPT invariant 7 t 7 rt (with R onsite) when 
the topological flux changes by /. In Ref. j5| it argued that 
when 7 ^, = — 1 , rjftj, = 1 (where R = c r in the notation of 
the Kagome model), there may be gapless edge modes pro¬ 
tected by a non-trivial vison PSG R~ 1 T~ 1 RT = —1 (as¬ 
sociated with topological superconductivity in the fermionic 
mean-field ansatz). Presumably the vison PSG can also be 
computed from the change in the ID SPT invariant 7 rt ( 
77- 1 for the vison) when threading a vison. 


VI. INTRINSIC TOPOLOGICAL ORDER: DETECTING 
THE TOPOLOGICAL FLUX 

It is important to have a method for measuring the topolog¬ 
ical flux a of an MES independent of the SET order; in the 
finite cylinder, we must detect whether a vison v threads the 
cylinder, and for a torus / infinite cylinder we must distinguish 
between all of 1, b, v, f. It has previously been shown that the 
topological S and T matrices can be calculated from the MES 
on both the torus and infinite cylinder, which can then be used 
to label the MES.|,33il In practice this is not so simple as all 
four-sectors must be found on a finite circumference cylin¬ 
der, which is frustrated by finite size effects which lead to a 
sizable splitting of the topological degeneracy. However, for a 
Z 2 spin-liquid we find there is a simpler procedure to uniquely 
label the MES individually. 

Note that an S' = 1/2 model must have finite topological- 
flux per unit cell li28ll . so the topological flux through an entan¬ 
glement cut depends on the location of the cut. So to simplify 
the discussion, we restrict to even circumferences L y £ 2Z, 
which contain an even number of S = 1/2 within each ring of 
the cylinder, and always consider ’vertical’ entanglement cuts 
which lie between ID super-sites under the 2D to ID dimen¬ 
sional reduction. 


A. {l,u} vs {6,/} 

The first distinction detects the fractional S = 1/2 spin 
carried by the spinons b and /. Under the dimensional re¬ 
duction, an even circumference cylinder is an SO(3) invariant 
integer-spin chain, which has a Z 2 ID SPT classification as¬ 
sociated with the emergence of two-fold degenerate S' = 1/2 
edge states (protected either by time-reversal or SO(3)). The 
simplest non-trivial example is the S = 1 AKLT state. In 
2D, when topological flux b or / terminates at the edge of the 
cylinder it also produces a spinon excitation near the edge car¬ 
rying S = 1/2. If the system is SO(3) symmetric, this emer¬ 
gent edge spin carries a two-fold degeneracy. Hence under the 
dimension reduction, the b, f sectors are non-trivial ID SPT 
states under SO(3), while 1, v are trivial. 


We conclude that the 1, v sectors will have an entangle¬ 
ment spectrum that transforms under integer representations 
of SO(3), and hence will have singlets in the entanglement 
spectrum, while the entanglement spectrum of the b, f sectors 
will transform under half-integral representations of SO(3), 
leading to a minimum two-fold entanglement degeneracy. 


B. 1 vs v 

In an SO(3) symmetric spin liquid, threading a vison v 
through the system is topologically equivalent to threading 
27 t - flux with respect to S z spin-rotations. To distinguish 
between the 1 and v sectors we can detect the change in mo¬ 
mentum induced by the flux threading. Viewing a cylinder of 
length L x and circumference L y as a spin chain on a periodic 
ring of length L y , the unit cell of the chain contains integral / 
half-integral spin when L x is even / odd. In the half-integral 
case, threading flux with respect to S z spin rotations is known 
to increase the y-momentum by e I7r = — 1. [.54] In the integral 
case, the threading flux will not change the momentum. 

We conclude that threading a vison through the system will 
increment the 7 -momentum by 7 r when L x is odd, and by 0 
if L x is even. This increment is simply the y-momentum per 
unit length r] X y V introduced in Sec. fill B I Hence for the vac¬ 
uum, rj xy = 1, while for the vison, r\ v xy = —1. This result is 
straightforward to check in any parton construction. 

We know r/ xy = —y/ y , since they differ by vison insertion, 
but which of the two carries r] xy = — 1 will in fact depend on 
the PSG relation (T x T y ) a = p a xy {T y T x ) a . 

C. b vs / 

A fermionic any on has topological spin 9f = —1, which 
we can use to distinguish between the b and / spinons. We 
will show that the fermion’s topological spin is encoded in an 
additional 2x degeneracy in the entanglement spectrum (ES); 
combined with the 2 x spin degeneracy, the fermionic spinon 
ES if 4x degenerate. Intuitively, the ES of the fermionic sector 
should have anti-periodic boundary conditions, meaning that 
the momenta k are quantized as k £ |^(Z + i). With ei¬ 
ther time-reversal or reflection symmetry, the momenta k and 
— k will be degenerate, so there is a 2x degeneracy. We will 
show this 2x degeneracy arises from non-trivial ID SPT order 
under a combination of Zr rotational symmetry and either 
reflection or time-reversal. 


1. Review of momentum polarization 

Momentum polarization is a procedure to detect topological 
spin using a translation T y that rotates a cylinder of circum¬ 
ference L y . 13(v 371 To review, each left Schmidt state |a) of 
the Schmidt decomposition {e ~ Ea , | a )} can be assigned def¬ 
inite momentum , meaning that T y \ a) = e lka \ a). In a 
convention in which k a £ /^-Z, the momentum polarization 
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is 

) Ly 

= T a e- {a ~ ir]H)L2 y + 0(e~ L y/t), 

(82) 

j- a = e 2*i(h a -c/2l) ' ( 83 ) 

Here a is non-universal real constant related to the area-law 
fluctuation of momentum of the cut; t]h is the ‘Hall-viscosity’ 
(which isn’t quantized on a lattice), and T a = e 2 ' Kl ( ha ~ c / 24: ' > is 
the desired entry of the modular T-matrix. 

The momentum polarization is not necessarily a ID SPT 
invariant, since there are no ID SPTs associated with an onsite 
1,L y symmetry. As a consequence, generally the momentum 
polarization only becomes quantized in the L y —> oo limit, 
and a scaling analysis is required. 

2. A ID SPT invariant for detecting fermionic topological flux 

However, in the presence of a mirror reflection y aa —y or 
time reversal, we can prove that the momentum polarization 
is a Z 2 ID SPT invariant that detects whether the topological 
flux is bosonic or fermionic. Consider a cylinder of even cir¬ 
cumference L y with a mirror reflection R y : (x,y) -A (x,—y) 
that acts as an onsite symmetry in the ID picture (the result 
for time reversal is analogous). When acting on the entangle¬ 
ment spectrum, a 7 r-rotation (translation by L y /2) may anti¬ 
commute with the inversion R y : 

U Ru (U Tv ) L '' / 2 U£ = (-l) F (UT y ) L > /2 , F = 0,1 (84) 

giving a ID Z 2 invariant F = 0,1. To show F is a ID SPT 
invariant, note the symmetry group generated by T y ,R y is 
G = Zi y x Z 2 (for even L y ). The cohomology classifica¬ 
tion is 

R 2 [Z Ly *Z 2 ,[/(1)]=Z 2 . (85) 

The relation F of Eq. ( 184b is a gauge invariant, so must label 
these two possibilities. 

To relate the SPT invariant F to the anti-periodic bound¬ 
ary conditions of the entanglement spectrum, suppose we 
redefine the phase of Up to ensure the expected relation 
Uf^U T yU R y = (f/r H ) _1 . With this gauge choice Eq. d84b 

requires (UT y ) Ly = (—1) F . For F = 1, the diagonal basis 
Ut v = e lka requires k a £ |^(Z + \). 

We note that the most general commutation relation is 

Ur v Ut v U^ = 1f U^ (86) 

with 7 p y ^ 2 = (— 1 ) F . 

3. Quantization of momentum polarization by the ID SPT 
invariant 

We now prove that the momentum polarization is quantized 
to be (—1) F ; this confirms the interpretation that F detects the 
topological spin 9f — — 1 of the fermionic spinon. 


Since Up y commutes with the entanglement spectrum 
e~ Ea , we have 

A = Tr(e- E "U Ty ) = Tr (e~ E -U^U Ty U Ry ) (87) 

= Tr(e" F “7 F (7 7 : 1 ) = Tr(e“ F “C/^) = 7 F A* (88) 

Using this relation, for even L y the momentum polarization is 

X Ly = X Ly/2 { lF X*) Ly/2 = fp y,2 \X\ Ly = {-l) F \X\ Ly . 

(89) 

So long at |A| 0, the momentum polarization is equivalent 

to the SPT invariant. The SPT invariant is even more robust 
since it is well defined even when |A| =0. 


VII. DETECTING SET ORDER ON THE KAGOME 
LATTICE USING CYLINDER-DMRG 

We now propose a procedure to determine the Kagome 
PSGs which is practically adapted to the constraints of cylin¬ 
der DMRG. The results of this analysis will be reported in a 
subsequent work.!59ll 

Following the notation of the earlier analysis, 0 within 
the fermionic parton construction we must determine the 
five invariants {r] a , r)c 6 , VaCe , V aT, VC 6 t, m2 = VC 6 r]aCe }> 
combined with a possible v x / y £ Z classification of any 
symmetry-protected gapless edge states. JJ] 


A. Finite DMRG 

We first consider a technique for finite length cylinders. By 
finding the SO(3) invariant ground state after adding or re¬ 
moving an extra spin at each edge of the geometry, DMRG 
studies can reliably obtain two topological sectors that differ 
by threading a spinon through the bulk. At the circumferences 
that can currently be well converged (such as YC8 and XC8), 
the other two topological sectors are not generally observed. 
Even without determining whether these two sectors are l/v , 
6 //, the ratio of R x , R y , //, quantum numbers before and after 
adding the extra sites (i.e. a spinon) will determine j , 

and referring to Table HU three of the PSG invariants. Since the 
symmetries must be edge-exchanging, on the YC type cylin¬ 
ders we obtain (/p/ 1 and hexagon centered ■ On the 

XC type cylinders, we obtain and hexagon centered 

Qjk/7), thg latter serving as a double check on the YC data. 
There is a simple algorithm for measuring space-group quan¬ 
tum numbers in finite DMRG.|59tl 

To determine ?yi 2 one must first determined whether the 
topological sector of the spinon is b/f. To distinguish b/f, 
we check if there is a 4-fold degeneracy in the entanglement 
spectrum as predicted by Eq. ( 184b . which would imply the 
sector is /. Knowing the sector b/f, we know the correct 
boundary condition for the parton ansatz, and following the 
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techniques of Sec. II V El we can predict the site-centered inver¬ 
sion Q Ia (A, b/ /). The details depend on the cylinder used, 
but the result always reveals 7712 . 

Once we know the topological sector of the sample, we 
can further check these results by comparing the absolute 
quantum numbers under R Xl R yi Ih/s with the computation 
of Sec. II V El there are many different cases depending on the 
sample. 

B. Infinite DMRG 

An infinitely long cylinder can be studied using iDMRG, 
which has certain numerical advantages due to the absence of 
edge effects and the reduced computational costs. iDMRG 
also reliably finds two topological sectors which differ by a 
spinon. Here we discuss only even-circumference cylinders. 

As discussed, the Z 2 ID SPT invariants for SO(3) / T de¬ 
termine which state has the spinon. The fermionic ID SPT 
invariant distinguishes between b/ f. The momentum per unit 
length r/zy is trivial to compute in iDMRG, so we distinguish 
between sectors 1/v using 1 = rjl y = — rj v xy . The momen¬ 
tum per unit length of the spinon sector determines either 
7712 = rjl y (if the sector is /) or —7712 = r] b xy (if the sector 
is b). 

To measure the reflection PSGs, one can in principle 
detect the ID SPT reflection invariants using established 
methods. 0 This can be a bit unwieldy in 2D DMRG, as 
the ordering of the DMRG ‘snake’ breaks the reflection sym¬ 
metries. A simpler procedure is to generate a finite cylinder 
wavefunction by projecting the left / right regions of the in¬ 
finite cylinder onto reflection related classical product states, 
leaving behind a finite segment of spins. [0] Regardless of 
projection used, the resulting state is a reflection symmetric 
finite cylinder wavefunction. One can then measure the space- 
group quantum numbers of the resulting finite cylinder wave- 
function in order to determine Ih 

C. Determining RTR~ 1 T~ 1 

As discussed in Sec. IV Cl the remaining invariants are re¬ 
lated to the onsite ID SPT invariants 7 rt for R = R x ,R y , 
which can be measured on YC and XC type cylinders respec¬ 
tively using known methods for detecting ID SPTs. 101 The 
most obvious signature is the 4-fold degeneracy in the ES re¬ 
quired to realize the projective relations Ut Of = 7 r = -1 
and {UrY^U^UtU^JJ^ 1 ) = "/t^rt = —1- The fermionic 
PSGs are related to the relative ID SPT order 7 (7) between 
the 1 and / sectors via 

= VcrTVcr (90) 

(f) (f) 

Tt Tr't = IvTlCeTriariaCs (91) 

( f) 

Of course 7 ^ = — 1 , as it is a spinon. 

If the DMRG obtains sectors which differ by /, we are 
done. If DMRG obtains sectors which differ by b, the analysis 


depends on whether the vison has a non-trivial RTR~ 1 T~ 1 
PSG. If the vison PSG is trivial, there are no gapless edge 
states and the boson will have the same RTR~ 1 T~ 1 PSG as 
the fermions, so the bosonic relative SPT order 7 ^ again re¬ 
covers the fermionic PSG invariants. If the vison PSG is non¬ 
trivial, there are gapless edge states, and tjaTVcr = — 1 for the 
R x edge and rjaTVCeTVvVcrCe = -1 for the R y edge regard¬ 
less. 


VIII. CONCLUSIONS 

We argued that the many body symmetry quantum numbers 
are a robust global property ideally suited to detecting distinc¬ 
tions between SETs, and calculated these for several SETs 
described by slave particle mean field theories with different 
projective symmetry groups. 

More generally, the SETs may be diagnosed from the ID 
SPT invariants in the cylinder geometry. These invariants for 
different topological sectors (which are labeled by the quasi¬ 
particles) combine together in a way that reflects the fusion 
rules. In contrast, combining PSGs for a pair of quasiparticles 
to predict the PSG for the fusion product needs to be carefully 
considered in the case of internal symmetries. 

The knowledgeable reader may be puzzled by this corre¬ 
spondence between PSGs and ID SPTs. The latter is de¬ 
termined by projective representations modulo a phase, or 
technically 'H 2 {G : (7(1)), where G is the symmetry group, 
while the former is a projective representation modulo Z 2 , 
'H 2 {G : Z 2 ), for Z 2 topological order. That is, we represent 
group elements by matrices, whose product satisfies the group 
multiplication, up to either a U(l) phase or just an overall sign 
(Z 2 ) of the matrices . This is because the physically observ¬ 
able quantities are made by combining two identical quasipar¬ 
ticles and so we can only change the overall sign of the ma¬ 
trices. For example, if X is a Z 2 symmetry, then ‘half charge’ 
of a quasiparticle corresponds to the PSG X 2 = —1. How¬ 
ever, this is not a ID SPT invariant. How is this discrepancy 
reconciled? 

The key observation is if inversion, I, or a equally a reflec¬ 
tion is present, these act like anti-unitary symmetries when 
restricted to the Schmidt states on one side of a bipartition. 
Thus, while I 2 = — 1 may again seem to be a PSG relation, 
regarding it as an antiunitary symmetry turns it into a projec¬ 
tive representation even with (7(1) phase factors and hence a 
ID SPT invariant UHHU]. Similarly, for a global Z 2 sym¬ 
metry, X, while by itself X 2 = —1 does not produce a ID 
SPT invariant, when combined with the effective antiunitary 
inversion symmetry ■ - - x x G — is a ID SPT invariant and the 
denominator is the fractional charge we are interested in. This 
type of reasoning has been repeatedly used in this work. 

Our procedure is expected to be complete for Z 2 liquids, but 
for more complicated topological orders, such as say Z 3 topo¬ 
logical order, there is a Z 3 invariant associated with a C 3 rota¬ 
tion symmetry. Taking into account certain subtleties, global 
C 3 quantum numbers detect this topological invariant, but the 
simplest cylinder dimensional reduction will not work. Exten¬ 
sions are left to future work. 









15 


IX. ACKNOWLEDGEMENTS 

We thank Steve White, David Huse and Zhenyue Zhu for 
a stimulating collaboration on related topics. We also thank 


Yang Qi and Liang Fu for alerting us to their recent preprint 
before submission, which overlaps with the present worklroltl. 
AV was supported by NSF-DMR 1206728 and the Templeton 
Foundation. MZ was supported by NSF DMR-1206515 and 
the David & Lucile Packard Foundation. 


[1] X.-G. Wen, Phys. Rev. B 65, 165113 (2002) 

[2] A. Kitaev, Annals of Physics 321, 2 (2006) january Special 
Issue. 

[3] A. Mesaros and Y. Ran, Phys. Rev. B 87, 155115 (2013) 

[4] A. M. Essin and M. Hermele, Phys. Rev. B 87. 104406 (2013) 

[5] Y.-M. Lu, G. Y. Cho, and A. Vishwanath, (2014), 
arXiv: 1403.0575 

[6] L.-Y. Hung and Y. Wan, Phys. Rev. B 87, 195103 (2013) 

[7] A. Vishwanath and T. Senthil, Phys. Rev. X 3, 011016 (2013) 

[8] L. Fidkowski, N. Lindner, and A. Kitaev, “unpublished,”. 

[9] F. Burnell, X. Chen. L. Fidkowski, and A. Vishwanath, (2014), 
arXiv: 1403.6491 

[10] M. Barkeshli, M. Bonderson, M. Cheng, and Z. Wang, (2014), 
arXiv: 1410.4540 

[11] H. Song and M. Hermele, (2014), arXiv: 1409.1504 

[12] X. Chen, Z.-C. Gu, and X.-G. Wen, Phys. Rev. B 83, 035107 

( 2011 ). 

[13] F. Pollmann, A. M. Turner, E. Berg, and M. Oshikawa, Phys. 
Rev. B 81, 064439 (2010). 

[14] N. Schuch. D. Perez-Garcia, and I. Cirac, Phys. Rev. B 84, 
165139 (2011). 

[15] L. Fidkowski and A. Kitaev, Phys. Rev. B 83, 075103 (2011) 

[16] X. Chen, Z.-C. Gu, Z.-X. Liu, and X.-G. Wen, Science 338, 
1604 (2012). 

[17] M. Levin and Z.-C. Gu, Phys. Rev. B 86, 115109 (2012) 

[18] Y.-M. Lu and A. Vishwanath, Phys. Rev. B 86, 125119 (2012) 

[19] P. A. Lee, Science 321, 1306 (2008) 

[20] L. Balents, Nature 464, 199 (2010) 

[21] J. S. Helton, K. Matan, M. P. Shores, E. A. Nytko, B. M. 
Bartlett, Y. Yoshida, Y. Takano, A. Suslov, Y. Qiu, J.-H. Chung, 
D. G. Nocera, and Y. S. Lee, Phys. Rev. Lett. 98, 107204 
(2007). 

[22] Y. Ran, M. Hermele, P. A. Lee, and X.-G. Wen, Phys. Rev. Lett. 
98, 117205 (2007). 

[23] Y. Huh, L. Fritz, and S. Sachdev, Phys. Rev. B 81. 144432 
( 2010 ). 

[24] O. Cepas, C. M. Fong, P. W. Leung, and C. Lhuillier, Phys. 
Rev. B 78, 140405 (2008). 

[25] S. Yan.D. A. Huse, and S. R. White, Science 332, 1173 (2011) 

[26] H.-C. Jiang, Z. Wang, and L. Balents, Nature Physics 8, 902 

( 2012 ). 

[27] S. Depenbrock, I. P. McCulloch, and U. Schollwock, Phys. 
Rev. Lett. 109, 067201 (2012). 

[28] M. Zaletel and A. Vishwanath, (2014), arXiv: 1410.2894 

[29] M. Levin and X.-G. Wen, Phys. Rev. Lett. 96, 110405 (2006) 

[30] A. Kitaev and J. Preskill, Phys. Rev. Lett. 96, 110404 (2006) 

[31] H. Li and F. D. M. Haldane, Phys. Rev. Lett. 101, 010504 
(2008). 

[32] E. Keski-vakkuri and X.-G. Wen, International Journal of Mod¬ 
ern Physics B 07, 4227 (1993). 


[33] Y. Zhang, T. Grover, A. Turner, M. Oshikawa, and A. Vish¬ 
wanath, Phys. Rev. B 85, 235151 (2012) 

[34] M. Heidar and X. G. Wen, (2014), arXiv:1401.0518 

[35] Y. Zhang, T. Grover, and A. Vishwanath, (2014), 
arXiv:1412.0677 

[36] M. P. Zaletel, R. S. K. Mong, and F. Pollmann, Phys. Rev. Lett. 
110, 236801 (2013). 

[37] H.-H. Tu, Y. Zhang, and X.-L. Qi, Phys. Rev. B 88, 195412 
(2013). 

[38] J. Haah, (2014), arXiv: 1407.2926 

[39] L. Cincio and G. Vidal, Phys. Rev. Lett. 110, 067208 (2013) 

[40] W. Zhu, S. Gong, F. D. M. Haldane, and D. N. Sheng, Phys. 
Rev. Lett. 112, 096803 (2014). 

[41] S.-S. Gong, W. Zhu, and D. N. Sheng, Sci. Rep. 4 (2014), 
10.1038/srep06317. 

[42] B. Bauer, L. Cincio, B. P. Keller, M. Dolfi, G. Vidal, S. Trebst, 
and A. W. W. Ludwig, (2014), arXiv:1401.3017 

[43] M. P. Zaletel, Phys. Rev. B 90, 235113 (2014) 

[44] F. Wang and A. Vishwanath, Phys. Rev. B 74. 174423 (2006) 

[45] Y. Huh" L. Fritz, and S. Sachdev, Phys. Rev. B 81. 144432 
( 2010 ). 

[46] S. Sachdev, Phys. Rev. B 45, 12377 (1992) 

[47] Y.-M. Lu, Y. Ran, and P. A. Lee, Phys. Rev. B 83, 224413 

( 2011 ). 

[48] F. Yang and H. Yao, Phys. Rev. Lett. 109, 147209 (2012) 

[49] C.-Y. Huang, X. Chen, and F. Pollmann, Phys. Rev. B 90, 
045142 (2014). 

[50] L. Wang, A. Essin, M. Hermele, and O. Motrunich, (2014), 
arXiv: 1409.7013 

[51] A. M. Turner, Y. Zhang, and A. Vishwanath, Phys. Rev. B 82, 
241102 (2010). 

[52] C. N. Varney, K. Sun, M. Rigol, and V. Galitski, Phys. Rev. B 
84, 241105 (2011). 

[53] S. Jiang, A. Mesaros, and Y. Ran, Phys. Rev. X 4, 031040 
(2014)." 

[54] E. Lieb, T. Schultz, and D. Mattis, Annals of Physics 16, 407 
(1961). 

[55] M. Oshikawa, Phys. Rev. Lett. 84, 1535 (2000) 

[56] Hastings, M. B., Europhysics Letters 70, 824 (2005) 

[57] Note that only the torsor structure of the MES is required here, 
which is why the procedure still applies to odd-circumference 
cylinders even though our notation is imprecise. 

[58] D. Perez-Garcia, M. M. Wolf, M. Sanz, F. Verstraete, and J. I. 
Cirac, Phys. Rev. Lett. 100, 167202 (2008) 

[59] Z. Zhenyue, M. P. Zaletel, Y.-M. Lu, D. Huse, A. Vishwanath, 
and S. White, “unpublished,”. 

[60] F. Pollmann and A. M. Turner, Phys. Rev. B 86, 125441 (2012) 

[61] Y.QiandL. Fu, (2015), arXiv: 1501.0009 



